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r ~7- Abstract 

We clarify important physics issues related to the recently established new mass 
limit for magnetized white dwarfs which is significantly super-Chandrasekhar. The 

Mh, issues include, justification of high magnetic field and the corresponding formation 

of stable white dwarfs, shape of the white dwarfs under consideration, flux freezing, 
variation of magnetic field and related currents therein. We also attempt to address 

j^ ■ the observational connection of such highly magnetized white dwarfs. 



> 1 Introduction 

Recently, Mukhopadhyay and his collaborators have proposed that highly magnetized 
white dwarfs could have significantly super-Chandrasekhar masses [H [21 [31 HI [5] . Notably, 
Das & Mukhopadhyay [I] have shown that such white dwarfs have a limiting mass of 
2.58M . In order to arrive at this result, the authors have made certain assumptions 
for the convenience of calculation. These are, the choice of constant magnetic field in 
the central region of white dwarfs, which decreases with density near the surface of the 
white dwarf, employing the idea proposed by previous authors in the context of strongly 
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magnetized compact objects [6]; the choice of spherical symmetry for strongly magnetized 
white dwarfs; separate central and surface flux conservation in magnetized white dwarfs 
as they evolve (e.g., due to accretion), etcetera. 

Now, in the present report, we plan to justify the above assumptions in a greater detail, 
in the light of various related physics issues. This is particularly important, as over the 
last few months several friends, colleagues and critics have raised certain questions about 
our proposal of significantly super-Chandrasekhar white dwarfs and a new mass limit of 
white dwarfs. This report is by and large a suitable compilation of our responses to about 
a dozen referee reports, which we have successfully tackled during this series of works 
starting with the one that appeared in 2012 [I]. 



2 Equipartition magnetic field for super-Chandrasekhar 
white dwarfs 

In this section, we would like to clarify that for the super-Chandrasekhar white dwarfs 
considered in our works, the equipartition magnetic field (B) is no longer same as that 
for weakly magnetized white dwarfs, but much higher. 

One can estimate the equipartition field of a star by recalling the scalar virial theorem 
[7]. By equating the magnetic and the gravitational potential energies of a spherical star 
of mass M and radius R, i.e., (4nR 3 /3)(B^ -/Sir) ~ GM 2 /R, one arrives at the maximum 
possible equipartition magnetic field -B equ i ~2x 10 s (M / M Q ) (R/ R Q )~ 2 G, where M® and 
R® are the mass and radius of Sun respectively and G is Newton's gravitation constant. 
For a typical white dwarf having radius ~ 10 4 km and mass 1AM®, this turns out to be 
~ 10 12 G. However, for a super-Chandrasekhar white dwarf having, for example, mass 
2.58M and radius ~ 70 km (see, e.g., Figure 1 of [5]), B cqui is ~ 5 x 10 16 G > 10 12 G. 
It is to be noted, however, that this is but an average value and hence is less than the 
value of the field considered to model the central region of the above white dwarf, and 
larger than the value at the surface. In fact, if the central region is defined to be the 
region of radius Ri n t = R/10 which demarcates a zone up to which the central magnetic 
field of the white dwarf remains (almost) constant [5], then the above argument yields 
-Bequi ~5x 10 18 G in this zone. This is very well consistent with the central magnetic field 
Bi n t (= 8.8 x 10 17 G) used to model the white dwarf considered above. Hence, although 
we are considering super-critical (central) magnetic fields (B > B c = 4.414 x 10 13 G), 
which cause Landau quantization to significantly alter the equation of state, we are not 
actually considering super-equipartition fields. Note that the value of -B equ i itself has 
increased for a super-Chandrasekhar white dwarf, since with the increase in the magnetic 
field strength, the white dwarfs start becoming more compact in size (see, e.g., Figures 
4(b), (c) and §IV.B of [2]). We would like to point out here that, indeed, Ostriker and 
Hartwick [5] constructed models of white dwarfs with B > 10 12 G at the center but with 
a much smaller field at the surface. Thus, the concept of high interior magnetic fields in 
white dwarfs, although difficult to verify observationally, are not implausible. 

3 Justification of constant central magnetic field 

In our works, we have modeled the central region of the white dwarf to have a strong 
constant magnetic field (several orders of magnitude larger than the surface field). This 
can be justified as follows. 

From Maxwell's equations in a steady state, we have 

An 
V x B = -j, (1) 

where j is the current density. Now, in the presence of a very strong central magnetic 
field, B ^> B c (as is the case for a significantly super-Chandrasekhar white dwarf), the 
electron degenerate matter will not only be Landau quantized but also the electrons will 
occupy mostly the ground Landau level. Thus the electrons are not expected to move and 
hence practically there will be no current, i.e., j = 0. Therefore, from the above equation, 
we obtain that there is no spatial variation of B, i.e., the magnetic field is constant. Since 



it is this central magnetic field which is primarily responsible for the mass to exceed the 
Chandrasekhar limit (see, e.g., §IV.C of [2J), one can consider the hydrostatic equilibrium 
condition in order to solve for the structure of the magnetized white dwarfs, at least in 
the central region, for obtaining super-Chandrasekhar white dwarfs. 

4 Stable equilibrium in presence of strong magnetic 
field 

One might wonder whether magnetic energy could really dominate over internal energy 
without magnetic forces significantly altering the equilibrium. Since the central region of 
the white dwarf in our works is modeled to have a constant magnetic field, the magnetic 
forces vanish therein. Nevertheless, the field must drop off by the time the surface is 
reached and hence the field gradients (and therefore magnetic forces) could exist primarily 
in the lower density outer regions, which are, however, less important for the super- 
Chandrasekhar mass calculation (see, e.g., Figure 5 and §IV.C of [2]). We demonstrate in 
this section that gravitational forces in the low density regions would still dominate over 
the magnetic forces making it possible for the white dwarf to attain a stable equilibrium. 

Let us consider a super-Chandrasekhar white dwarf having a mass 2.39M (which 
is less than the mass limit of 2.58M ), radius 9 x 10 7 cm and a central magnetic field 
Bint — 4.414 x 10 15 G (this is the white dwarf corresponding to the maximum mass point on 
the dotted line in Figure 2 of [4]). Figure 1(a) shows the variation of density as a function of 
radius inside the above white dwarf having a central density p c = 4. 14 x 10 9 gm/cc. Figures 
1(b) and 1(c) show the variation of matter pressure and mass as functions of radius. If we 
consider a radius in the low density region of this white dwarf R' — 7 x 10 7 cm, then from 
the figures we see that the density at this radius is p(R') = 1.46 x 10 9 gm/cc, the mass 
within this radius is M(R') = 1.73M© and the matter pressure is P(R') = 8.83 x 10 26 
ergs/cc. Adopting the description of the variation of magnetic field with density in a 
compact star given by Equation 5 of Bandyopadhyay et al. [6] , we show a plausible profile 
of the field in Figure 1(d) below, using suitable values for the parameters (3 = 0.001 and 
7 = 6 (note, it is not the same as the polytropic index) and the values of central and 
surface magnetic fields. We note that the profile assures almost constant B in the central 
region of the white dwarf and then B drops off near the surface. Since it is the central field 
which is responsible for contributing significant mass to the white dwarf, such a profile 
would produce very similar results to those we have already obtained. From Figure 1(d) 
we obtain that the magnetic field at p(R') is B(R') = 1.6 x 10 14 G. 

Now, in the presence of a varying magnetic field the white dwarf must obey the 
condition of magnetostatic equilibrium: 
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where, r is the radial distance from the center of the white dwarf, P and p are the pressure 
and density of the electron degenerate matter respectively. We now present order of 
magnitude estimates of the various terms of Equation (J2J) in the low density region of the 
white dwarf, where the gradient of magnetic field is non-zero. The orders of gravitational, 
pressure and magnetic forces per unit mass, at R ', are GM(R')/R' 2 , P(R')/(p(R')R') and 
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Figure 1: (a) Density, (b) pressure and (c) mass, as functions of radius, within a white 
dwarf of mass 2.39M Q and central magnetic field 4.414 x 10 15 G. (d) Variation of magnetic 
field with density inside the same white dwarf, po is a normalization density with a value 
of 8 x 10 8 gm/cc. R is in units of cm, p in gm/cc, P in ergs/cc, M in M and B in G. 



B (R') 2 J (SttR' 'p(R')) respectively. Hence at R' the gravitational, pressure and magnetic 
forces per unit mass are of the order 4.71 x 10 10 dyne/gm, 8.64 x 10 9 dyne/gm and 
9.96 x 10 9 dyne/gm respectively. Thus we see that the gravitational force still dominates 
over the magnetic force (and combined rest of the forces) in the low density region and 
hence it is possible to construct a stable equilibrium solution. Note that, when exact 
gradients are calculated in reality, gravitational force will exactly balance the combined 
forces due to matter pressure and magnetic field. 

5 Degenerate matter density stronger than magnetic 
density 

One can also make stability arguments by computing the magnetic density (ps = B 2 /87rc 2 ) 
and the matter density (p) inside the white dwarf shown in Figure 1 and then their effects 
on the gravitational energy density. Now, the mass of the white dwarf depends directly 
on the matter density and the gravitational energy density is in turn proportional to the 
mass. Hence, higher the matter density, higher is the gravitational energy density. The 
magnetic density and the matter density at the center of the white dwarf of mass 2.39M Q 
are 8.6 x 10 8 gm/cc and 4.14 x 10 9 gm/cc respectively. Thus the matter density is about an 
order of magnitude larger than the magnetic density. The magnetic field remains constant 
in the central region (and hence so does ps) up to a matter density ~ 3 x 10 9 gm/cc (see 
Figure 1(d)), but still we note that matter density is greater than ps = 8.6 x 10 8 gm/cc. 
As the density decreases further, the magnetic field starts falling off near the surface of the 
white dwarf. For example, the magnetic densities pb = 2.1 x 10 8 gm/cc, 3.5 x 10 6 gm/cc 
and 7.8 x 10 3 gm/cc correspond to the matter densities p = 2.4 x 10 9 gm/cc, 1.6 x 10 9 
gm/cc and 9.6 x 10 8 gm/cc, respectively. Thus we see that, ps decreases much faster than 
p as the surface of the white dwarf is approached and hence a stable equilibrium can be 
maintained throughout the white dwarf and ps can be ignored to compute the mass of 
white dwarfs. Note that, if in the central region the contribution from ps is included in 
order to compute the mass of the white dwarfs, then this may slightly change the mass 
estimate and probably even increase it. 

6 Stability against magnetic buoyancy 

A question that may arise is the issue of buoyancy related instabilities within the strongly 
magnetized super-Chandrasekhar white dwarfs. If a blob of gas is displaced adiabatically 
from an inner radius to an outer radius in the central region of the white dwarf, then the 
condition for stability is given by: 



j I K i j, 



i/r 2 

"' = {¥/') > *■ (3 > 

where p' is the matter density inside the blob when it is at the outer radius, pi, T\ and 
K\ are the matter density, adiabatic index and a dimensional constant related to the 
adiabatic equation of state respectively of the surrounding medium at the inner radius 
and P2, T 2 and K 2 are the same quantities at the outer radius. Note that, since the 
magnetic field is constant in the central region of the white dwarf, the magnetic pressure 
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Figure 2: Variation of the difference in densities inside and outside a displaced blob of 
gas as a function of its surrounding density. Density is in units of gm/cc. 



terms do not appear in the above relation. The above criterion means that the blob at 
the outer radius must be heavier than the surrounding (i.e., p' — p2 > 0) so that it is 
pulled inwards by gravity to its original position, making the system stable. If the blob 
is lighter than its surrounding (i.e., p' — p<i < 0) then the buoyant force displaces the blob 
further outwards, making the system unstable. 

We apply this argument to a white dwarf having central density ~ 1.4 x 10 9 gm/cc; 
the equation of state of the underlying degenerate gas being described by the parameters 
listed in Table II of [2] . In Figure 2 we show the variation of p' —p2 with p 2 inside this white 
dwarf in a few discrete radial zones. We clearly see that p' — p 2 is everywhere positive, 
establishing that the white dwarf is stable against buoyancy related instabilities. Away 
from the central region, the magnetic field decreases towards the surface (as discussed 
previously), so the magnetic pressure at the outer radius will be less than that at the 
inner radius. Taking this difference into account will only increase the density of the blob, 
making it even more heavier than the surrounding, further strengthening the argument 
for a stable system. 



7 Choice of spherical symmetry 

Order of magnitude estimates obtained in §4 argue that the overall gravitational energy 
in the white dwarfs under consideration dominates over the combined total energy that 
arises from matter and magnetic field. This justifies our choice of spherical symmetry. 



However, if the magnetic field is very strong, then the pressure may become anisotropic 
[2]. Following Bocquet et al. [TO], the combined fluid-magnetic medium then may develop 
a tension at such a high magnetic field. Bocquet et al. [TO] observed that as a result of 
this magnetic tension, the star displays a pinch across the symmetry axis and assumes 
a flattened shape. It is possible that a similar effect could occur in highly magnetized 
white dwarfs. In that case, the magnetic tension would deform the magnetized white dwarf 
along the direction of the magnetic field and possibly turn it into a kind of oblate spheroid. 
Consequently, the spherically symmetric hydrostatic equilibrium equation would have to 
be modified due to the presence of the strong magnetic field. As stated by Paulucci et 
al. [11], one has to generalize the Tolman-Oppenheimer-Volkoff equation to take into 
account the anisotropy in the pressure, but to the best of our knowledge this has not 
been achieved so far, as is also mentioned in the same paper itself. We propose to carry 
out such calculations for highly magnetized white dwarfs in a future work. However, in 
the previous work [2] , we have initiated the understanding of the effect of deviation from 
spherical symmetry due to the magnetic field, if any, in the Appendix. Interestingly we 
observe that, flattening effects due to magnetic field renders super-Chandrasekhar white 
dwarfs even at relatively weaker magnetic fields which are more probable in nature. 

Important point to note, however, is that, for the strongly magnetized white dwarf 
having the new limiting mass [3], the above geometrical factors do not matter, since the 
corresponding radius is negligible (as is also the case for the Chandrasekhar limit). 

8 Difficulty behind direct observation of strongly mag- 
netized super-Chandrasekhar white dwarfs 

Recent observations of several peculiar over-luminous type la supernovae - SN 2006gz, SN 
2007if, SN 2009dc, SN 2003fg - seem to suggest super-Chandrasekhar-mass white dwarfs 
with masses up to 2.1 — 2.8M as their most likely progenitors [T21 HH]- The strongly 
magnetized super-Chandrasekhar white dwarfs proposed by Das & Mukhopadhyay [21 H] 
and Das, Mukhopadhyay & Rao [5] can explain these supernova explosions, which is their 
indirect evidence. However, there has been no direct detection of such white dwarfs so 
far. We try to explore the reason behind this in this section. 

Ordinary sub- Chandrasekhar white dwarfs are located in the lower-left portion of the 
Hertzsprung-Russell (H-R) diagram (plot showing the relationship between the luminosi- 
ties of stars and their effective temperatures), indicating that they are low luminosity hot 
objects. If the surface of the white dwarf is considered to radiate like a black body, then 
its luminosity is given by L = AuFPcjbT^^ where gb is the Stefan-Boltzmann constant 
and T e ff the effective surface temperature of the white dwarf. Since L oc R 2 , ordinary 
white dwarfs, having radius ~ 10 4 km, are bright enough to feature in the H-R diagram. 
Neutron stars on the other hand, having radius 10 km, are too small (and hence too 
dim) and would lie far below the white dwarfs in the H-R diagram, and hence do not 
not appear in it altogether. Note that, neutron stars are mostly detected as pulsars, due 
to the emissions of regular pulses. The strongly magnetized super-Chandrasekhar white 
dwarfs in our works, having radius < 100 km, are more closer to neutron stars than con- 
ventional white dwarfs. Thus, such (isolated) white dwarfs would also not appear in the 
H-R diagram, making them virtually impossible to be detected directly based on current 
observational techniques. 



9 Separate flux conservation in the central and sur- 
face regions 

Das, Mukhopadhyay & Rao [5] have discussed a possible evolutionary scenario from weakly 
magnetized, sub-Chandrasekhar, accreting white dwarfs to strongly magnetized, super- 
Chandrasekhar white dwarfs, which eventually explode to give rise to the peculiar type la 
supernovae. Irrespective of how the interior magnetic field of the white dwarf is modeled, 
the only constraint on the above evolution is the magnetic flux conservation between the 
initial and final white dwarfs, which controls/relates the magnetic field profiles of the 
white dwarfs, if one of them (say the initial one) is known. Additionally, the authors have 
imposed that the magnetic flux is conserved separately in the central and surface regions 
of the white dwarfs. This is justified as follows. 

In §2 of [5], the authors have defined a radius Ri n t which demarcates a zone up to 
which the central magnetic field of the white dwarf remains (almost) constant which is 
larger than that at the surface. This is quite a plausible assumption, since the original 
star which collapses to form a white dwarf might have a very large interior magnetic field 
compared to that observed on its surface. For example, Parker [T3] estimated that a 
magnetic field strength as high as ~ 10 7 G can not be ruled out for the Sun's interior. 
Gough & Mclntyre [15] also discussed the inevitability of a strong magnetic field in the 
Sun's interior and its influence to the extent of a transitional layer, called the tachocline. 
Thus, it is quite likely that the values of the central and surface magnetic fluxes (in the 
initial star and hence in the white dwarf) and the underlying electric currents might be 
significantly different. Therefore, the size of the central (and then the surface) region of 
a white dwarf is controlled by the size of the central region of the original star. Since it is 
likely that there is a gradient of magnetic field from the center to the surface of the star, 
the radial dependence of B could be of the form B oc 1/r 1 *, where ( is positive. However, 
the above field profile will blow up at r = which is unfeasible and hence B is expected to 
be very large and to remain constant in the central region. This restricts the electrons to 
move around in the central region of the final highly magnetized white dwarfs, as only the 
lower Landau level(s) is/are available for them, justifying practically no current therein. 
Nevertheless, as B decreases towards the surface, higher Landau levels become available 
and hence the electrons can produce current, maintaining a gradient of the magnetic field 
up to the surface. 

All the above arguments justify the assumption of the separate conservation of the 
surface and central magnetic fluxes, which also suggest the presence of a core-envelope 
boundary. 

10 Correlation between the observed peculiar type 
la supernovae rate and the occurrence of highly 
magnetized super-Chandrasekhar white dwarfs 

Identification of the progenitors of type la supernovae is still an ongoing research area. In 
[5] we have proposed that if the progenitors are single degenerate white dwarfs accreting 
matter from a normal companion, then the peculiar supernovae originate from super- 
Chandrasekhar white dwarfs with extremely high magnetic fields. 
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The observed rate of over-luminous, peculiar type la supernovae is about 1-2% of the 
normal type la [16J. This is consistent with the fact that about 25% of all cataclysmic 
variables (CVs) are magnetic CVs and we expect a small fraction (say 10%) of them 
to be super-magnetic. These high magnetic CVs (2.5% of all CVs based on the above 
assumption) would behave like normal magnetic CVs during most of their life time, till 
their mass remains sub-Chandrasekhar. However, they would be seen as peculiar objects 
only during a small fraction of their life time, when the surface magnetic field becomes 
very high (> 10 11 G) and the mass is super-Chandrasekhar (> 1.44M ), which should 
happen because of accretion. We have discussed this issue with constant and varying 
accretion rates in §3 of [5]. 

For a constant accretion rate of 2.5 x 10 _9 M Q yr _1 , which leads to a white dwarf having 
final mass 2.33M (shown by the dashed line in Figure 3(a) of [5]), the white dwarf 
spends about 0.4 fraction of its lifetime with super-Chandrasekhar mass. For varying 
accretion rate, modeled by Equation 4 of [5] with a = 2 (solid line in Figure 3(b) of 
[5]), the white dwarf spends about 0.05 fraction of its lifetime with super-Chandrasekhar 
mass. Now, typically observed stably accreting intermediate polars have accretion rates 
~ lO _9 M yr _1 [T7]. However, an accretion rate > lO _7 M yr _1 is required for the stable 
burning of the accreted hydrogen and helium in (single degenerate) type la supernovae 
progenitors [TS], which is indeed attained by the above model with varying accretion rate, 
starting with an accretion rate similar to that of intermediate polars. Thus, the varying 
accretion rate is a more likely scenario, which gives us the estimate that about 0.12 % of 
CVs would be peculiar high magnetic CVs. Since they are accreting at a high rate, they 
will be visible as bright X-ray objects, rather than white dwarfs in optical data bases like 
SDSS. 

11 Summary 

We have reestablished the subtle physics issues related to highly magnetized super- 
Chandrasekhar white dwarfs. We have shown that such white dwarfs are indeed possible, 
provided certain obvious physics constraints are satisfied, explained in this report. This 
further strengthens the argument for the existence of such super-Chandrasekhar white 
dwarfs. 
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